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ABSTRACT 


The  Goncept  of  equivalent  surface  reaction  is  developed  for  gas-phase 
recornbination  at  the  stagnation  region  Of  blunt  bodies.  This  concept 
is  based  on  the  fact  that  the  chemical  state  of  a  nonequilibrium  highly 
cooled  boundary  layer  is  largely  determined  by  the  recombination  that 
occurs  near  the  wall.  An  equation  based  on  this  concept  is  shown  to 
predict  the  heat  transfer  to  a  nohcatalytic  surface  to  within  a  few 
percent  of  the  more  accurate  existing  results.  The  equation  is  gener¬ 
alized  to  apply  to  the  case  in  which  surface  catalytic  recombinatioh 
occurs  simultaneously  with  gas -phase  reaction. 
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INTRODUCTION 


Ttie  chemical  state  of  nonequilibrium  boundary  layers  about  hypersonic 
vehicles  is  of  considerable  interest.  The  interest  stems  mainly  from  the  fact 
that  the  heat  transfer  to  the  body,  and  the  observables  in  the  wake  and  the 
trails,  are  greatly  affected  by  the  chemical  state  of  the  boundary  layers. 

The  problem  of  nOnequilibrium  chemically  recombining  boundary  layers 

is  a  formidable  one  and  exact  solutions  are  available  only  for  a  few  special 

cases.  The  solutions  that  exist  have  been  obtained  by  assuming  that  either 

1  -  6* 

the  gas  phase  or  the  surface  recombination  occurs  alone.  At  this  stage 

of  development,  an  approximate  solution  will  be  useful  if  it  takes  reasonable 
account  of  the  simultaneous  effects  of  gas  phase  and  surface  recombination. 

7 

It  has  been  shown  that  the  chemical  state  of  the  boundary  layer  is 
controlled  mainly  by  the  recombination  which  takes  place  near  the  wall  when 
the  wall  is  highly  cooled.  This  observation  suggests  representing  the  gas 
phase  reaction  by  an  equivalent  surface  reaction  with  all  state  variables 
specified  by  surface  conditions.  The  order  of  reaction  and  the  corresponding 
rate  constant  of  this  equivalent  surface  reaction  will  be  determined  for  a 
noncatalytic  wall.  The  equivalent  surface  reaction  concept  will  then  be 
extended  to  the  case  in  which  the  surface  has  an  arbitrary  catalycity. 

The  basic  concept  of  emphasizing  the  region  of  maximum  chemical 

activity  in  determining  the  chemical  state  of  the  boundary  layer  was  first  used 

8 

by  Rae  for  the  boundary  layer  in  which  dissociative  reactions  predominate. 

ANALYSIS 

Fay  and  Riddell^  have  obtained  an  exact  solution  of  the  nonequilibrium 

stagnation  point  boundary  layer  with  the  chemical  state  of  the  gas  in  the 

boundary  layer  ranging  from  frozen  to  equilibrium.  The  computations  were 

for  a  noncatalytic  wall  at  a  temperature  of  300°K  and  for  a  selected  flight 

2 

condition.  Chung  and  Anderson  have  calculated  by  an  integral  method  the 

Numbers  indicate  References  at  the  end  of  paper. 


nonequilibriiim  effect  on  heat  transfer  to  blunt  cones  with  noncatalytic  walls 

at  temperatures  of  1000  and  15O0°K  for  two  flight  conditions.  Goodwin  and 
9 

Chung  have  successfully  obtained  additional  Splutibns  and  have  correlated 
the  results  for  various  flight  GOnditions  and  a  wall  temperature  of  1500°K  by 
a  Single  recombination  rate  parameter. 

7 

Subsequently,  Inger  has  shown  that  the  difference  in  wall  temperature 

1  2  9 

between  the  two  theoretical  solutions  ’  ’  may  be  included  in  the  correlation 

by  noting  the  controlling  influence  of  the  wall  temperature  upon  the  chemical 
state  of  the  boundary  layer.  It  appears  that,  with  the  resulting  GOmprehensive 
correlation,  the  essential  effects  of  nOnequilibrium  reaction  on  chemical 
state  and  heat  transfer  have  been  taken  into  account,  at  least  to  the  correct 
order  Of  magnitude. 

The  order  and  rate  constant  for  the  proposed  equivalent  surface  reaction 
will  be  determined  first  for  noncatalytic  walls.  Since  the  gas  phase  reaction 
near  the  Surface  of  a  highly  cooled  hypersonic  vehicle  is  predominantly  that  of 
recombination  requiring  three  body  collisions,  a  modified  second  order 
reaction  appears  to  be  a  reasonable  choice  for  the  equivalent  surface  reaction. 
We  let 


/8  m 
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w\  8  y 


1  = 


m 


K 


1  + 


(1) 


where  is  the  equivalent  surface  reaction  constant  for  the  gas  phase  reaction 
The  constant  is  to  be  determined  subsequently.  In  effect,  we  are  assuming 
that  the  nonequiiibrium  gas  phase  reaction  over  a  noncatalytic  wall  is  equiv¬ 
alent  to  a  frozen  boundary  layer  with  the  surface  reaction  specified  by  Eq.  (1). 

.  3 

An  exact  solution  of  the  diffusion  equation  at  the  stagnation  point 
gives  the  following  equation  for  a  frozen  boundary  layer 

(|E)  =0.47(1  (2) 
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Solution  of  Eqs.  (1)  and  (2)  for  gives* 
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For  the  Liewis  number  not  much  different  from  one,  it  can  be  readily 


q  -  qf 


=  1 
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Equations  (4)  and  (6)  give  at  the  stagnation  point 


-11^  =  1  -  /HZ 

'lE  -  ^f  V  1  +  K 

''  g 


For  the  stagnation  region: 

u  =  Px 
e 

r  =  X 


(3) 


(4) 


(5) 

shown  that 

(6) 


(7) 
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Solutions  of  nonequilibrium  boundary  layers  exist  for  axi symmetric  bodies 

only.  The  following  correlations,  therefore,  will  be  limited  to  the  stagnation 

9 

regions  of  axisymmetric  bodies.  The  correlation  Of  Goodwin  and  Ghung 

7 

including  the  surface  temperature  correlation  Of  Inger  gives 
(q  -  qjj/fqg  ^  qj)  in  terms  of  a  Single  parameter: 


r 


g 


It  is  desired  to  find  such  that  the  heat  transfer  ratio  (q  -  q£)/(qj,  -  q£), 
obtained  by  Eq.  (7)  will  match  the  result  given  by  Goodwin  and  Chung  in  terms 
of  r^.  It  is  found  that  the  heat  transfer  ratios  match  quite  satisfactorily  if 

k  =  21  r  (9) 

g  g 


Figure  1  shows  the  comparison  between  the  heat  transfer  ratio  obtained  from 
Ref.  9  and  that  from  Eqs.  (7)  and  (9).  It  is  seen  that  the  present  equivalent 
surface  reaction  theory  represented  by  Eqs.  (1),  (7),  and  (9)  predicts  the  heat 
transfer  ratio  satisfactorily  for  noncatalytic  wails. 

Let  us  now  consider  surfaces  of  arbitrary  catalycity,  represented  by  a 
specific  rate  constant  K^,  in  the  presence  Of  nonequilibrium  gas  phase 
reaction.  The  diffusion  rate  of  atoms  to  the  surface  must  then  be  equal  to  the 
total  rate  of  recombination  of  atoms  by  both  the  gas-phase  recprnbinatipn  and 
the  catalytic  surface  recombination.  Equation  (1)  is  generalized  for  this 
case  as 


D  (^) 
r  w\oy/. 


C  —  +  {p  K  )m 

g  1  +  m  "^w  w  w 


The  solution  of  Eqs.  (2),  (10),  and  (6)  for  (q  -  -  q£)  with  the 
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Md  of  Eqs.  {5.)  and  (9)  gives 


q  -  q>  f  (r  +  4(1  +  r  +  2ir  -  r 
^  _  ,  L  w _ _ w  g  J _ w 

■tj,  ‘  4j  ■  ■  EW  r,  ■»  2iy  ~ 


where  for  the  stagnation  region  Of  axisymrfietric  bodies 


r 
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(p  K 

Tw  w _ 

(0.47)[2(p^p^)^BG]l/^ 


(11) 


(12) 


Equation  (11)  is  a  closed  form  solution  which  gives  the  heat  transfer  ratio 
(q  -  q^)/(q^  -  q£)  for  arbitrary- rate  simultaneous  gas-phase  and  surface 
reactions . 


The  terms  r  and  r  .  defined  by  Eqs.  (8)  and  (12),  are  essentiaUy 
g  w  ’ 

Damkholer  numbers  representing  the  ratio  of  characteristic  diffusion  time 
to  characteristic  chemical  reaction  time  for  the  gas -phase  and  surface 
reactions,  respectively. 

DISCUSSION 

The  accuracy  of  the  solution,  Eq.  (11),  cannot  be  checked  for  the 

arbitrary  combinations  of  and  because  the  general  exact  solutions  are 

not  available.  It  can  be  established  that  the  equation  gives  correct  results 

at  the  two  limiting  cases  of  either  F^  =  G  or  F^  =  0.  It  is  sufficiently  accurate 

when  F  =0  because  F  was  obtained  by  matching  the  results  with  the  solutions 
w  g 

of  Ref.  9.  It  is  also  correct  when  F^  =  0  because  the  solution  (11)  essentially 
becomes  the  solution  of  Ref.  3.  Though  there  is  no  direct  check  in  the 
intermediate  region,  it  is  felt  that  the  physical  reasoning  leading  to  the 
equivalent  surface  reaction  concept  of  Eq.  (10)  is  sound.  The  excellent 
correlation  shown  for  the  gas  phase  reaction  with  noncatalytic  surface 
appears  to  support  the  accuracy  of  the  results  presented. 
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The  variation  of  (q  ~  qj')/(q^  -  q^)  with  respect  to  an4  F^,  as  given 
by  Eq.  (li)  is  seen  in  Fig.  1.  The  general  behavior  is  similar  to  that  for 
the  Couette  flow  Obtained  in  Ref.  10. 

It  is  noted  here  that  the  same  problem  considered  here  is  also  being 
studied  by  Inger*  from  a  different  approach.  In  his  analysis  the  production 
term,  or  the  source  term,  in  the  formal  diffusion  equation  is  approximated 
by  a  simpler  function  requiring  only  a  correct  matching  near  the  surface, 
it  is  therefore  based  on  a  concept  similar  to  that  used  herein  in  that  the 
chemical  state  near  the  surface  is  assumed  to  be  of  controlling  influence. 

It  is  still,  however,  an  approximate  analysis  and  the  accuracies  of  any  of 
these  approximate  analyses  can  be  checked  only  when  exact  solutions  become 
available . 

GONCI.USION:S 

The  concept  of  equivalent  surface  reaction  has  been  developed  for  gas- 
phaSe  recombination  at  the  stagnation  region  of  blunt  bodies.  This  concept  is 
based  on  the  fact  that  the  chemical  state  of  a  nonequilibrium  highly  cooled 
boundary  layer  is  largely  determined  by  the  recombination  that  occurs  near 
the  wall-  An  equation  based  on  this  concept  was  shown  to  predict  the  heat 
transfer  to  a  noncatalytic  surface  to  within  a  few  percent  of  the  more  accurate 
existing  results.  The  equation  was  generalized  to  apply  to  the  case  in  which 
surface  catalytic  recombination  occurs  simultaneously  with  gas-phase  reaction. 
The  present  solution,  which  is  in  a  simple  closed  form,  should  be  useful  in 
estimating  heat  transfer  although  the  accuracy  of  the  solution  cannot  be 
checked  for  the  general  case,  due  to  a  lack  of  exact  solutions.  The  atom 
concentration  profile  across  the  boundary  layer  can  be  obtained  from  the 
value  m^  obtained  herein,  since  in  a  frozen  boundary  layer  the  profile  is 
determined  when  m  is  known  (for  instance,  see  Ref.  3)  and  the  present 
theory  is,  in  essence,  based  on  the  frozen  boundary  layer  theory  with  modi¬ 
fication  only  in  the  boundary  condition  at  the  surface.  The  profile  thus 


*T6  be  published. 
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obtained  should  a.pproximate  the  true  profile  with  the  same  degree  of  accuracy 
as  the  heat  transfer  results.  This  profile,  in  turn,  defines  the  chemicai 
state  across  the  boundary  layer. 

The  equivalent  surface  reaction  concept  will  be  extended  tO  flow  fields 
beyottd  the  stagnation  point  of  a  blvint  body. 
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NOMENCiAtTJEE 

G  =  p|i  ratio  defined  by  Eq.  (3) 
e  =  total  atom  mass  fraction 

=  atom  mass  fraction  of  nitrogen  atoms 
D  =  binary  diffusion  coefficient 
=  frozen  total  enthalpy 
Ah°  =  heat  of  recombination 

=  equivalent  surface  reaction  constant  for  gas  phase  reaction 
=  parameter  defined  by  Eq.  (5) 

=  specific  catalytic  surface  recombination  constant 

20 

=  constant  portion  of  recombinatipn  coefficient  (=  1.  56  X  10  ) 

m  =  c/c 

e 

p  =  pressure  (atmosphere) 
q  =  heat  transfer  to  surface 

R  =  universal  gas  constant  (=  82.  06  cm^  atm/mole  °K) 
r  =  distance  from  the  axis  of  symmetry  to  the  surface 
s  =  function  defined  by  Eq.  (3) 

Sc  =  Schmidt  number,  p/pD 
T  =  absolute  temperature  {°K) 
u  =  streamwise  velocity 
X  =  streamwise  distance 
y  =  distance  normal  to  surface 
p  =  (du  /dx)  (in/sec  for  Eq.  8) 


( 


f ) 


Fg  -  DafiakiiSJJer  number  for  gas -phase  reGombination 

r  -  Damkhbler  number  fbr  surface  reebmbinatibn 
w 


'€  = 


j:0  for  two-dimenSiohal  body 
;  1  for  axisymnfietrie  body 
T1  -  funGtion  defined  by  Eij-  (3) 
p  -  viscosity 

V  =  kinematic  viscosity,  p/p 
p  -  density 


StXBSGRIiPfS 


E  -  equilibrium 
e  -  edge  of  boundary  layer 
f  =  frozen 

o  =  stagnation  point 
w  =?  wall 
oo  -  free  stream 
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